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Abstract. 

Although  topological  error-correcting  codes  offer  a  promising  paradigm  for  fault- 
tolerant  quantum  computation,  their  robustness  in  the  presence  of  leakage  to  non- 
computational  states  is  unclear.  Here  we  explore  the  signature  and  consequences 
of  leakage  errors  on  ancilla-assisted  Pauli  operator  measurement  in  superconducting 
devices.  We  consider  a  realistic  coupled-qutrit  model  and  simulate  the  repeated 
measurement  of  a  single  oz  operator.  Typically,  a  data-qubit  leakage  event  manifests 
itself  by  producing  a  “noisy”  ancilla  qubit  that  randomly  reads  |0)  or  |1)  from  cycle 
to  cycle.  Although  the  measurement  operation  is  compromised,  the  presence  of  the 
leakage  event  is  apparent  and  detectable.  However,  there  is  also  the  possibility  of  a  less 
typical  but  more  dangerous  type  of  leakage  event,  where  the  ancilla  becomes  paralyzed , 
rendering  it  oblivious  to  data-qubit  errors  for  many  consecutive  measurement  cycles 
and  compromising  the  fault-tolerance.  Certain  dynamical  phases  associated  with  the 
entangling  gate  determine  which  type  of  leakage  event  will  occur  in  practice.  Leakage 
errors  occur  in  most  qubit  realizations  and  our  model  and  results  are  relevant  for  many 
stabilizer-based  error  correction  protocols. 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 
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1.  Introduction 


Topological  quantum  error- correcting  codes,  such  as  surface  and  toric  codes,  are 
attracting  attention  because  of  their  high  error  thresholds  and  realistic  designs  that 
only  require  nearest- neighbor  interactions  [1,  2,  3,  4,  5,  6].  While  the  robustness  of 
standard  fault-tolerant  topological  codes  to  discrete  Pauli  errors  is  a  subject  of  active 
research  [6,  7,  8,  9],  the  effect  of  leakage  to  non-computational  states  still  remains  an  open 
question  and  is  relevant  for  most  quantum  computing  architectures.  Understanding  the 
effect  of  such  leakage  errors  is  important  for  superconducting  qubits  not  only  because 
higher  energy  states  |2)  ,  |3)  , . . .  are  present,  as  is  the  case  with  most  other  qubit 
realizations,  but  also  because  they  can  be  utilized  to  implement  two-qubit  entangling 
operations  such  as  the  Strauch  controlled-crz  (CZ)  gate  [10,  11].  One  way  to  suppress 
the  effects  of  leakage  is  to  adopt  the  topological  cluster-state  approach  [2,  3,  12],  where 
each  qubit  is  repeatedly  initialized,  operated  on  by  gates,  and  measured:  This  approach 
systematically  removes  leakage  errors  from  all  qubits  in  the  array,  at  the  cost  of  some 
additional  operations.  Another  approach  might  be  to  use  a  stabilizer-based  topological 
error-correcting  code  for  qudits ,  and  theoretical  progress  has  been  made  in  this  direction 
[13,  14,  15,  16,  17,  18,  19], 
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Figure  1.  Protocol  for  ancilla-assisted  az  measurement.  Here  “A”  is  the  ancilla  qubit 
and  “D”  the  data  qubit.  Each  cycle  (dashed  box)  consists  of  a  reset  of  the  ancilla  to 
the  |0)  state,  a  Hadamard  gate  H  on  the  ancilla,  a  CZ  gate,  and  another  Hadamard 
followed  by  ancilla  readout  in  the  diagonal  basis.  The  readout  result  is  recorded  and 
the  cycle  is  repeated  indefinitely.  The  data  qubit  never  gets  measured  or  reset. 


In  this  work,  however,  we  consider  the  standard  stabilizer-based  approach  for 
qubits — but  applied  to  three-level  qutrits — and  regard  any  population  transfer  to  the 
1 2)  state  of  the  ancilla  or  data  qutrit  as  a  potential  error.  We  investigate  the  origin  and 
signature  of  such  leakage  errors  for  an  ancilla-assisted  measurement  of  the  data  qutrit, 
identify  a  potentially  dangerous  regime  where  data  errors  are  invisible  to  the  ancilla, 
and  discuss  its  consequences  for  topological  error  correction. 

Figure  1  shows  the  circuit  for  our  protocol.  Let’s  review  how  this  works  in  the  ideal 
limit:  Initially,  the  data  qutrit  D  is  assumed  to  be  in  some  pure  qubit  (not  qutrit)  state 

IV’d)  =  a|0)  +  6|1) ,  (1) 

while  the  ancilla  A  is  initialized  to  |0).  We  perform  the  gate  operations  shown,  record 
the  measurement  outcome,  reset  the  ancilla  to  |0),  and  repeat  this  cycle  many  times. 
Throughout  this  work,  the  Hadamard  gate  (H)  is  assumed  to  be  ideal  and  to  act  as  the 
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identity  on  the  third  level  of  the  qutrit, 
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The  Hadamards  and  CZ  combine  to  produce  a  controlled-NOT  (CNOT)  gate  that  copies 
the  data  qubit  to  the  ancillla,  but  here  we  implement  this  CNOT  with  the  gates  shown 
in  Fig.  1  because  we  believe  that  the  CZ  gate  can  be  implemented  in  superconducting 
architectures  with  very  high  fidelity  [11],  For  an  initial  D  state  (1),  the  state  of  the 
system  after  the  second  H  gate  is,  in  the  |AD)  basis, 


a  |00)  +  b  1 11) . 


(3) 


Thus,  in  the  absence  of  any  errors,  the  readout  projects  the  data  qubit  into  the  observed 
eigenstate  of  the  ancilla.  And  once  the  data  qubit  is  projected  to  a  computational  basis 
state,  it  remains  there  forever. 

In  this  work  we  study  the  effects  of  intrinsic  gate  errors  and  decoherence  on 
this  process.  The  remainder  of  the  paper  is  organized  as  follows:  In  Sec.  2  we 
describe  our  physical  model  and  consider  ancilla- assisted  measurement  in  the  presence 
of  decoherence.  The  non- ideal  CZ  gate  is  discussed  in  Sec.  3.1.  Leakage  errors  and 
ancilla  paralysis  are  discussed  in  Sec.  3.2.  We  discuss  the  implications  of  our  results  for 
the  design  of  error-corrected  superconducting  quantum  computers  in  Sec.  4. 


2.  Coupled  qutrit  model 

In  this  section  we  describe  our  model,  and  for  a  warm-up,  show  how  the  ancilla-assisted 
measurement  protocol  works  with  ideal  gates,  but  in  the  presence  of  decoherence. 


2.1.  Model 


The  Hamiltonian  for  a  pair  of  capacitively  coupled  transmon  or  phase  qutrits  is  given 
by 
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where 


V  = 


0  -i  0 
i  0  —iy/ 2 
0  iy/2  0 


(5) 


Qutrit  1  is  the  ancilla  qutrit  and  qutrit  2  is  the  data  qutrit.  In  (5)  we  have 
assumed  harmonic  qutrit  eigenfunctions.  The  time-dependence  of  the  Hamiltonian  (4) 
is  embedded  in  the  qubit  frequencies  €\  and  62;  the  Hadamard  gates  are  implemented 
with  microwaves  via  terms  not  shown  in  (4).  For  the  CZ  gate  protocol,  we  assume  the 
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Figure  2.  (Color  online)  Energies  of  various  levels,  in  the  |AD)  basis,  as  a  function 
of  ei/27 r.  Here  £2/277  =  6  GHz,  the  coupling  strength  is  <7/277  =  25 MHz,  and 
771/277  =  772/277  =  200  MHz. 


frequency  of  the  data  qubit  to  be  fixed  at  6  GHz,  while  the  ancilla’s  frequency  is  varied. 
The  anharmonicities  77^ / 27t  are  assumed  to  be  equal,  frequency-independent,  and  fixed 
at  200  MHz.  Figure  2  shows  the  energies  of  several  relevant  eigenstates  as  a  function 
of  e  1 ,  with  62/2-77  =  6  GHz  and  coupling  strength  g/2-K  =  25  MHz.  Note  that  the  only 
anticrossing  at  6i  =  e2  +  771  (ei/27r  =  6.2  GHz  in  Fig.  2)  is  between  the  1 11)  and  1 20) 
channels;  we  use  this  anticrossing  for  our  CZ  gate  [10,  11]. 

The  CZ  gate,  both  ideal  and  non-ideal,  is  parameterized  in  this  work  via  its 
generator.  A  generator  of  any  unitary  matrix  U  is  defined  as  a  Hermitian  matrix  S 
such  that  U  =  elS .  For  a  two-qutrit  system,  the  generator  of  the  ideal  CZ  gate  is  a 
Hermitian  matrix  S,  whose  matrix  representation  in  the  basis 

AD)  =  {  00) ,  01)  02) ,  |10) ,  |11) ,  1 12) ,  20) ,  |21) ,  |22)  } 


(6) 
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Note  that  within  the  computational  subspace,  elS  acts  as  a  standard  CZ  gate,  while 
non-computational  basis  states  acquire  phases  e1^ .  We  emphasize  that  any  extension  of 
an  ideal  CZ  gate  to  qutrits  is  dependent  on  the  assumed  model  and  gate  protocol.  For 
the  Strauch  CZ  gate,  auxiliary  crz  rotations  on  the  ancilla  and  data  qubits  nullify  the 
phases  acquired  by  the  |01)  and  1 10)  channels.  Since  we  use  the  anticrossing  between 
1 11)  and  1 20) ,  they  acquire  a  phase  of  angle  7 r.  We  assume  that  the  gate  is  in  the 
adiabatic  regime,  and  that  the  parameters  £*  are  dynamical  phases,  which  can  then  be 
expressed  as 
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00  0 

Here  tgate  is  the  operation  time  for  the  CZ  gate  (including  auxiliary  2  rotations),  and 
Eij  is  the  energy  of  eigenstate  \ij),  shown  in  Fig.  2.  To  keep  our  analysis  general  we 
do  not  assume  specific  values  for  the  £*.  They  depend  on  the  details  of  the  CZ  gate 
implementation  but  remain  fixed  throughout  a  given  experiment  or  simulation  (unless 
one  changes  tgate  or  the  pulse  shape).  As  we  will  explain  below,  the  difference 


d  =  £2  -  £1 


^gate 


€1  dt 


(9) 


Jo 

determines  if  the  ancilla  becomes  paralyzed  during  a  leakage  event.  Note  that  9  can  be 
varied  during  an  experiment  by  changing  the  gate  time. 
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2.2.  Ancilla-assisted  measurement  with  decoherence 

As  shown  in  Fig.  1,  each  measurement  cycle  consists  of  ancilla  initialization,  three  gate 
operations,  and  ancilla  readout.  Assuming  ideal  gates,  the  data  qutrit  after  the  first 
cycle  is  projected  to  a  computational  |0)  or  |1)  state  depending  on  the  observed  state 
of  the  ancilla  [recall  (3)].  In  the  absence  of  any  errors,  the  measurement  outcome  of  the 
ancilla  remains  unaltered  thereafter.  However,  the  situation  is  different  in  the  presence 
of  decoherence. 

In  order  to  model  the  effects  of  decoherence  on  the  measurement  outcomes  of  the 
ancilla,  we  assume  that  the  readout  and  reset  operations  are  instantaneous,  while  the 
Hadamard  and  CZ  gates  take  10  and  25  ns  respectively.  We  also  assume  that  amplitude 
damping  is  the  only  source  of  decoherence,  in  which  case  the  single-qutrit  Kraus  matrices 
can  be  written  as 
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For  an  operation  of  time  duration  At, 

Xm  =  l-e~mAt/T\  (11) 


cycle 


Figure  3.  (Color  online)  Simulated  repeated  readout  of  the  ancilla  qutrit  in  the 
presence  of  amplitude  damping.  Single  peaks,  upward  or  downward,  indicate  errors  on 
the  ancilla.  Data  errors  result  in  steps;  an  example  is  shown  near  cycle  1000.  In  this 
simulation  we  assume  T±  =  40  /xs,  T2  =  2Ti,  and  tcyc\e  =  45  ns. 

We  simulate  the  ancilla-assisted  measurement  protocol  for  an  ideal  CZ  gate  but  in 
the  presence  of  decoherence,  for  40,000  consecutive  cycles,  and  Fig.  3  shows  a  typical 
outcome.  The  duration  tcycie  of  each  complete  measurement  cycle  is  45  ns  (one  CZ  gate 
plus  two  Hadamards).  Initially,  the  data  qutrit  is  in  state  1 1) ,  and  a  single  downward 
peak  denotes  an  error  on  the  ancilla.  Near  the  1000th  cycle  the  data  qutrit  relaxes  to  |0) 
due  to  decoherence,  and  once  in  the  ground  state  it  stays  there  forever.  The  remaining 
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upward  peaks  are  caused  by  decoherence  on  the  ancilla  qutrit.  Since  the  ancilla  gets 
reset  at  the  end  of  every  cycle,  such  errors  are  manifested  as  single  peaks.  Note  that  if 
the  initial  state  of  the  two-qutrit  system  is  inside  the  computational  subspace,  it  does 
not  leak  to  non-co mputational  states  and  therefore  Fig.  3  is  insensitive  to  the  values  of 
the 


3.  Non-ideal  CZ  gate 

In  this  section,  we  first  discuss  how  a  non-ideal  CZ  gate  is  parameterized  and  then 
investigate  its  action  on  the  ancilla- assisted  qubit  measurement. 


3.1.  Parameterization  of  the  non-ideal  CZ  gate 


Let  us  first  give  a  brief  review  of  the  dominant  intrinsic  error  mechanisms  that  are 
relevant  for  the  Strauch  CZ  gate  [11];  the  Hadamards  are  always  assumed  to  be  ideal 
[see  (2)].  The  CZ  gate  of  Strauch  et  al.  [10]  is  performed  by  using  the  anticrossing 
between  the  1 11)  and  1 20)  states  at  ei  =  e2  +  77 1.  Although  the  other  states  are  detuned 
from  each  other  at  this  anticrossing  point,  a  small  amount  of  nonadiabatic  population 
transfer  is  unavoidable,  and  these  nonadiabatic  excitations  dominate  the  intrinsic  gate 
errors.  These  errors  can  be  thought  of  as  producing  a  second  unitary  matrix  whose 
generator  S'  can  be  parameterized,  in  the  basis  (6),  as 
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•  (12) 


The  complete  non-ideal  CZ  gate  is 

TT  _  i(S+S') 

Ucz  —  e  , 


(13) 


where  S  is  the  generator  (7)  of  the  ideal  CZ  gate.  The  parameters  Xi  and  Ci  in  (12) 
are  small,  while  the  angles  0*  take  arbitrary  values  between  0  and  2n.  Ci  and  C2 
parameterize  the  errors  occurring  during  pre  and  post  az  rotations,  and  C3  and  C4  denote 
the  controlled-phase  error  for  the  1 1 1)  and  |20)  channels.  In  our  simulations  we  assume 
Xi  =  Ci  —  10'2  for  all  i  —  1, ...  ,4.  Because  population  transfer  probability  scales  with 
| Xi  ,  our  choice  of  parameters  bounds  the  intrinsic  gate  errors  to  about  10-4. 


3.2.  Leakage  events  and  ancilla  paralysis 

The  CZ  gate  (13)  produces,  on  any  1 1 1)  input  component,  a  small  amplitude  of  1 02) 
(the  amount  determined  by  X2)  and  1 20)  (determined  by  x3)-  A  1 20)  component  either 
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results  in  the  possibility  of  an  ancilla  readout  of  |2) — if  the  readout  protocol  distinguishes 
|1)  and  |2) — or  the  possiblity  of  an  isolated  ancilla  error  if  it  does  not.  Neither  case 
compromises  fault-tolerance.  The  parameter  y2  is  responsible  for  data  qubit  leakage 
events.  By  a  leakage  event  we  mean  a  near-unity  population  of  the  data  \2)  state. 

The  principal  mechanism  producing  a  leakage  event  is  the  abrupt,  nonlinear 
transformation  on  the  data  qutrit  induced  by  the  ancilla  measurement.  We  denote 
these  transformations  by  T0,  Tl5  and  T2,  where  the  subscript  corresponds  to  the  ancilla 
readout  result.  Repeatedly  measuring  the  ancilla  applies  a  random  sequence  of  the  T 
maps  to  the  data  qutrit. 

For  the  model,  gate  implementation,  and  parameter  values  considered  in  this  work, 
the  map  T0  is  primarily  responsible  for  the  observed  leakage  events.  Although  the 
general  form  of  T0  is  quite  complex,  it  is  possible  to  construct  a  simple  special  case  of 
it  that  exhibits  the  essential  features.  To  do  this  we  choose  simplified  parameter  values 

6  =  7T, 

</>*  =  o, 

0  =0, 

Xs  =  0, 

X4  =  o,  (14) 

and  calculate  the  action  of  the  non-ideal  measurement  circuit  on  an  arbitrary  data  qutrit 
state 


IV’d)  =  a|0)  +  6|1)  +  c|2). 

We  find  (in  the  |AD)  basis)  that 


(15) 


»|00)  +  6|01)  +  e|02)  ->  |0>®  (|  +  ^  cos  y i  +  ^sinxi)|0) 

(b  a  .  b  c  .  \  |  .  ( c  b  . 

+  (j  cos  xi  -  ^  sm  Xi  -  ^  cos  X2  -  -  sm  X2J  1 1)  +  y -  +  -  sm  y2 

b 


cos  y2  1 2)  +|1) 


a  a 


L  V  2 


cos 


Xi  ~  g  sin  Xi)  1°)  +  (2  cosxi 


d  fo  c  \  /  c  b  c 

-  2  sin  Xi  +  2  cos  W  +  2  sin  X2J  | !)  -  y -  -  -  sin  y2  +  -  cos y2  )  |2)  .  (16) 

An  ancilla  readout  result  of  |0)  then  induces  the  map  T0  given  by 

a  +  a  cos  xi  +  b  sin  xi 


a  —$■  a  = 


b^b'  = 


b  cos  xi  —  a  sin  Xi  —  b  cos  y2  —  c  sin  y2 

: 


,  c  +  b  sin  X2  -  C  cos  X2 
c  — >  c  = - 


(17) 


Af  =  \a  +  acosxi  +  6  sin  %i  | 2  +  |6cosxi  —  a  sin  X;\  ~  b  cosy2  —  csiny2| 


+  |c  +  b sin X2  ~  c cos y2 


(18) 


where 
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Using  (17)  we  find  that  in  the  limit  Xi  —  0  and  y2  ~ t  0  the  data  qutrit  prepared  in  the 
|1)  state  transforms  as 

X0|l)  =  |2>.  (19) 


Our  simulations  confirm  that  the  dominant  mechanism  for  producing  a  leakage  event  is 
the  process  (19). 

Once  leaked,  the  data  qutrit  remains  in  the  |2)  state  (for  many  cycles)  until  it 
either  undergoes  a  nonadiabatic  “reverse-leakage”  transition  or  it  relaxes  back  to  the 
computational  subspace.  The  behaviour  of  the  ancilla  during  a  leakage  event  depends 
on  the  values  of  £1  and  £2  in  (7).  While  the  data  qubit  is  in  the  |2)  state,  the  two-qutrit 
system  is  restricted  to  the  subspace  spanned  by 


{ |02) , | 12) }, 


(20) 


because  the  1 22)  state  is  decoupled  and  remains  unoccupied.  In  this  subspace,  the  CZ 
gate  (13)  acts  as 


exp 


[  /  u 

0  ii 

H  0 

6  ). 

(21) 


and  therefore  performs  a  £  rotation  on  the  ancilla  by  an  angle  (9).  The  ffadamards  in 
Fig.  1  convert  this  to  an  x  rotation  [see  (9)] 

e-i{e,2)**  (22) 


acting  on  the  initial  ancilla  state  |0).  Therefore,  during  a  leakage  event,  while  the  data 
qubit  is  locked  in  the  |2)  state,  the  state  of  the  ancilla  after  every  cycle  is 

cos  1 10)  +  sin  1 1 1),  (23) 

and  upon  measurement  the  ancilla  qubit  reads  |0)  with  probability  cos2  (0/2). 

For  example,  if 

9  mod7r  =  ^,  (24) 

we  will  observe  random  ancilla  outcomes  with  equal  probabilities  for  observing  |0)  and 
|1).  This  type  of  leakage  event  is  simple  to  detect  (and  possibly  correct),  ffowever,  if 

6  mod  7T  =  0,  (25) 


then  the  ancilla  will  always  read  1 0) ,  cycle  after  cycle,  giving  no  indication  of  the  data 
error  and  thereby  compromising  fault-tolerance.  We  refer  to  this  dangerous  phenomena 
as  ancilla  paralysis. 

Figure  4  shows  the  readout  values  generated  from  the  sequential  measurements 
of  the  ancilla  qubit  for  different  choices  of  6 ,  including  all  error  process  contained  in 
the  non-ideal  CZ  gate  (13).  While  we  observe  random  oscillations  for  larger  values  of 
6,  no  such  signature  is  present  for  6  =  0.  In  order  to  quantify  the  paralysis  of  the 
ancilla  we  define  a  metric  W,  which  is  the  average  spacing — number  of  cycles — between 
consecutive  readouts  of  |1).  In  the  absence  of  decoherence,  we  can  estimate  it  [see  (23)] 


as 


W  =  csc2(0/2), 


(26) 
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Figure  4.  (Color  online)  Simulated  sequential  measurements  of  the  ancilla  qubit.  The 
readout  values  |0)  or  |1)  are  shown  as  a  function  of  measurement  cycle  number.  Red 
rectangles  signify  leakage  events,  where  the  data  \2)  state  probability  is  close  to  unity. 
Random  ancilla  oscillations  during  the  leakage  events  are  observed  except  when  9  i=s  0. 
Two  values  of  W  are  given  for  each  trace:  the  theoretical  value  from  (26)  and  a  value, 
shown  in  parentheses,  numerically  computed  from  the  simulation.  The  simulations 
assume  T\  =  40  fis,  T2  =  2Ti,  \i  =  Ci  =  10-2  for  alii  =  1, . . .  ,4,  and  random  values 
of  phase  angle  parameters  consistent  with  the  indicated  values  of  9. 


which  agrees  well  with  the  numerical  simulations. 

The  detectability  of  a  leakage  event  depends  on  whether  W  is  small  enough  to  be 
observed  in  the  presence  of  a  background  value  W*  resulting  from  decoherence  (and 
possibly  other  errors).  For  example,  in  the  simulations  of  Fig.  4,  which  have  Ti  =  40 
ps  and  T2  =  2Ti,  the  average  spacing  between  ancilla  |1)  peaks  away  from  the  leakage 
events  is  2381  cycles,  which  is  not  too  far  from  the  crude  theoretical  estimate 


W* 


2T\ 

/■cycle 


=  1778, 


(27) 


using  /cycle  =  45  ns.  The  estimate  in  (27)  can  be  derived  from  the  Pauli  twirling 
approximation  for  qubit  decoherence  [6] ,  which  predicts  ax  and  ay  errors  on  the  ancilla 
with  probability  px  =  pY  =  /Cycie/4Ti,  leading  to  a  total  bit-flip  probability  px  +  pY  of 
/cycle / 2T\ .  We  can  use  (26)  to  estimate  the  critical  value  of  9  separating  the  region  of 
dangerous  ancilla  paralysis  and  that  of  ordinary  leakage,  namely 


-l 


J 


Q*  —  2  esc  1  '/II'*  ss  2  esc 


(28) 
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which  is  9*  =  0.04  in  the  simulations  reported  here.  CZ  gates  with  9  mod7r  <  9*  are 
susceptible  to  undetectable  leakage  events. 

4.  Conclusion 

We  have  studied  the  basic  ancilla-assisted  measurement  circuit  in  the  presence  of  leakage 
errors,  and  identified  a  rare  but  potentially  dangerous  ancilla  paralysis  effect  that  could 
compromise  the  error-detecting  ability  of  a  stabilizer  measurement  operation.  Whether 
or  not  undetectable  paralysis  will  ocur  depends  on  the  difference  (9)  of  phase  angles 
produced  by  the  CZ  gate.  Although  fault-tolerance  is  compromised  with  either  type  of 
leakage  event,  the  ability  to  detect  such  an  event  might  allow  one  to  reset  the  affected 
qubit  to  recover  from  it.  We  note  that  the  value  of  9*  is  likely  to  be  larger  in  a  multi¬ 
qubit  Pauli  measurement  because  the  cycle  time  is  longer  (the  background  value  W*  is 
larger).  Our  results  suggest  that  leakage  be  addressed  either  at  the  hardware  level,  by 
periodically  removing  any  |2)  state  probability,  or  by  using  an  architecture  such  as  the 
2D  topological  cluster  code,  where  every  qubit  gets  measured  during  the  error-correction 
cycle.  In  addition,  it  is  of  course  advantageous  to  adjust  9  mod  tt  to  a  safe  value  near 
tt/2. 
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